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Creation of electron or exciton by external fields in a system with initially statistically independent 
unrelaxed vibrational modes leads to an initial condition term. The contribution of this term in the 
time convolution generalized master equation approach is studied in second order of the perturbation 
theory in the path integral formalism. The developed approach, applied for analysis of dynamics in 
the photosynthetic reaction center, exhibits the key role of the initial condition terms at the primary 
stage of electron transfer. 



PACS numbers: 87.15.ht, 05.60.Gg, 82.39.Jn 
I. INTRODUCTION 

The effect of environment on transport properties of 
quantum systems is a highly topical problem in atomic, 
nuclear and condensed matter physics. Markovian ap- 
proaches have been successfully used to study various 
phenomena in open systems, when the past memory of 
the system is neglected [l| . The advent in ultrafast laser- 
pulse technology [1] , quantum information processing Q , 
synthesis of new superheavy elements in cold and hot- 
fusion reactions |4|] are requiring a resolution of quantum 
dynamics, when a system is far from equilibrium. 

Analogous processes can be found in biological systems 
@ as well as in nanoscale devices Although there 
are wide structural and functional differences, the laws, 
that govern quantum solar energy conversion to chemi- 
cal energy or electricity in biological systems and semi- 
conductor solar cells, share many similarities. In these 
systems the conversion processes proceed from the cre- 
ation of electron-hole pairs (exitons) by a photon of light, 
followed by charge separation to produce the required 
high-energy product. The efficiency of solar cells may 
be increased due enhanced multiple exciton production 
in semiconductor quantum dots, which is essentially cre- 
ated extremely fast upon absorption of high-energy pho- 
tons [tI- On the other hand, the photosynthetic reaction 
center (RC) of bacteria provides an interesting system 
for studying a high-efficiency electron transfer in an or- 
ganized molecular complex. The RC is a special pigment- 
protein complex, that functions as a photochemical trap. 
In such systems, after excitation the electron transfer is 
so fast that there can exist unrelaxed vibrational modes 
in the primary stage of electron transfer. 

In a most studied case of a purple bacteria the RC 
is composed of three protein subunits called L, M and 
H [1, l^l- According to experimental facts, the pro- 
tein H does participate in the electron transfer (ET). 
All molecules (cofactors) involved in the ET are non- 
covalently bound to subunits L and M in two chains. 



Both chains of cofactors start at the bacteriochlorophyll 
dimer (P) which interacts with both subunits L and M. 
The dimer plays the role of the donor of an electron (a 
weakly bound exciton) at the photon absorption. Co- 
factors in the subunit L are accessory bacteriochloro- 
phyll (Bi), bacteriopheophytin (H^) and quinone (Ql). 
Identically, in the M subunit there are the accessory 
bacteriochlorophyll (Bm), bacteriopheophytin (H^) and 
quinone (Qa/). The cofactors serve as donor-acceptor 
pairs in the ET. The arrangement of cofactors shows the 
local twofold symmetry which is almost perfect with the 
respect to the dimer. The part of the L subunit involved 
in the ET can be superimposed onto the corresponding 
part of the M subunit by a rotation of almost exactly 180° 
(for more details on structural arrangement see jlO|). 

In spite of the structural symmetry of the two chains 
of cofactors, it appears that the RC is functionally highly 
asymmetric. In the primary charge transfer an electron is 
transferred from the photoexcited dimer P, the starting 
point for a series of electron transfer reactions across the 
membrane, to the cofactors on subunit L, to B^, H^, Q^, 
and Qm [HI [13 • On the other hand, the chain located on 
subunit M is inactive in the ET. The highly asymmetric 
functionality, however, can be decreased by amino acid 
mutations or cofactor modification [l^. If the direct ET 
between subunits L and M is not allowed and electron 
cannot escape from the system then it was shown that the 
different stochastic fluctuations in the energy of subunits 
and the interaction between subunits on these two ways 
may cause the transient asymmetric electron distribution 
at L and M branches during relaxation to the steady state 
[3]. However, due to fast electron transfer the memory 
effects should be important at the primary processes in 
the photosynthesis. The major goal of the present paper 
is to elucidate the effect of the initial condition terms on 
electron transfer in a system with initial conditions being 
far from equilibrium. 

The content of the paper is as follows. In Sec. II we 
outline the derivation of the generalized master equation 
adapted for physics of the RC. In Sec. Ill we use this 
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equation to analyse the contribution of the initial condi- 
tion terms on electron transfer in three-site model of the 
RC. A brief summary is presented in Sec. IV. In Appendix 
A there are details of the GME equation derivation. 

II. MASTER EQUATION FOR REACTION 
CENTER 

Formally, an exact generalized master equation (GME) 
which describes the electron transfer processes in systems 
with dissipation can be constructed by means of the pro- 
jection operator techniques To be specific let us 
consider a system in which an electron has N accessible 
sites embedded in a medium. Such a system is described 
by the Hamiltonian 

H = Ho + V, (1) 

where Hq = J2n=i W)[£n - «r„ -I- H^]{n\. Here, |n) is 
the electron state with energy e„. The parameter h/2Tn 
characterizes the electron lifetime at site \n). It may orig- 
inate, for example, from a nonradiative internal conver- 
sion or a recombination process. The term describes 
a medium (a solvent) consisting of harmonic oscillators 

= E + I'^^^lil^ - d„^A . (2) 

Here, dna is the equilibrium configuration of the ath os- 
cillator, which depends on the electronic state \n). The 
interaction V = J2n m=i ^nmln) {m\{n ^ m) couples dif- 
ferent sites. 

The GME formulation of an electron (an exciton) in- 
teracting with vibrational modes (phonons) starts from 
the Liouville equation for a density operator p{t) [l| 

i^^p{t) = i [Hp{t) - p{t)H+] EE Lp(t) . (3) 

At this stage the projector operator technique allows to 
avoid a knowledge of a complete information upon a sys- 
tem under consideration. The projector operator con- 
tracts the full information about the system to the rel- 
evant one. Our prime interest is the information about 
the electron localization and the irrelevant information is 
a particular vibrational state excited in the system. In 
virtue of this technique p^ . [l6| one obtains 

^Dp{t) = -iDLDp{t) - (4) 




DLeyi]y[-i{l - D)Lt]{1 - D)LDp{t - T)dT - 



-iDLeyi^[-i{l - D)Lt]{l - D)p{0) 

for the relevant part Dp{t) of the total density p(t). Here 
D = is an arbitrary linear projection operator which 
can be used in the form 

DA = Y^TT{\n){n\ A)pr,\n){n\. (5) 

n 



The total trace Tr = Tr^Trq is a product of traces of 
the electronic {Tr^} and the vibrational {Ttq) subsys- 
tems; pn = exp{—/3H^^)/TrQ[exp{—/3Hj^)] is a density 
operator for a vibrational subsystem, when an electron 
is localized at a site \n) (/? = 1//cbT). Often, the projec- 
tor operator is chosen in such a way that the initial state 
(1 — D)p{0) is disregarded This approximation is 
valid when electron and phonon subsystems are initially 
in equilibrium. Evidently, that at the primary stage of 
the electron transfer this term may influence the electron 
pathway in the RC. The questions arise about the time- 
scale of such influence in the system under consideration 
and how this influence would affect the ET. 

To calculate the initial condition terms we have to 
specify p(0). The photon absorption by the dimer re- 
sults in the transition of an electron from the ground 
to the excited state of the dimer (say, the excited state 
|1)). Before the excitation the electronic subsystem is in 
the thermal equilibrium with the vibrational subsystem 
which consists in vibrational modes of the dimer and the 
protein subunits. Due to the fast electron transfer to 
molecules located in L(M)-branches of the RC the time 
is too short to establish a thermal equilibrium between 
the vibrational subsystem and a new electronic state. On 
the other hand, the vibrational subsystem (the bath) is in 
the thermal equilibrium with the electronic ground state 
of dimer. Thus, we have 

p{0)^Po®(^pUO)\k){l\^ (6) 

where po — exp(— /^ffp )/TrQ[exp(— /?i7Q )] is a density 
operator for a vibrational subsystem, when an electron 
is at the ground state. Here p^;(0)=TrQ((fc|p(0)|Z)) is an 
electronic part of the density matrix. We suppose that 
an electron is initially localized on the first molecule: i) 
Pii(O) = Ij Pnn(O) = 0, n = 2, A^; u) non-diagonal initial 
density matrix elements are p^;(0) = 0. We also assume 
that after the excitation the electron transfer is so fast 
that the initial vibrational density is not affected. Sim- 
ilar assumptions for the construction of the initial state 
p(0) have been used for analysis of the energy transfer 
dynamics in a model of a donor-acceptor pair [isj . In 
principle, the initial conditions may be calculated within 
a scheme proposed to include system-bath correlations 
after the interaction with optical pulses [l^. However, 
this problem requires a dedicated study itself in order to 
distinguish different time scales and is beyond the scope 
of the present consideration. The main objective here is 
to gain insights into the role of the initial condition terms 
on the ET for a few typical cases. 
With the aid of Eq. Q one obtains 

i?p(0)=pi®|l)(l| , (7) 

which leads to the initial state 

(l-i?)p(0) = (po-pi)®|l)(l|. (8) 
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Evidently, if the condition po ~ pi is not fulfilled, one 
must take into account the initial state in Eq.Q. 
Substituting Eq.® in Eq.g]), one obtains the GME 



The ITs fulfil the general identity which transforms in 
the considered case to 



-ReWmn{t - T)Pra{T)]dT + 



(9) 



for site occupation probabilities 

Pn{t)=Tr{ \n){n\ p{t) ) = p„„(t) 



(10) 



The first term in the r.h.s. of Eq.® is associated with 
the probability for an electron to escape from the system 
via an additional channel. 

In the Born approximation, the memory function 
Wmn{i) can be expressed in the form 



W |2 



exp 



exp I 



ft2 

t I exp 



t X (11) 



E 



h J \ hiJa 

[(n„ + l)e-*'"°* + n„e''""* - (2n„ + 1)] 



Here ria — [6x^(11^^/ k^T) — 1] ^ is a thermal population 
of the ath vibrational mode and 



(12) 



is the reorganization energy of the ath vibrational mode, 
when an electron moves from state \m) to state \n). 

In second order of the perturbation theory, one obtains 
from Eg . ([4t the initial condition term (IT) in the form 
(see also [20, [HI) 

/„(<)== /«(i)+4^'(t) (13) 
where the first order term is 

/(i)(i) ^ -^^r( \n){n\Ce-^^^^\l~D)p{Q)) (14) 
and the second order term has the form 



/(2)(t) = - fdT^t,r). 
Jo 



(15) 



Here 



it,T) = Tr (|n)(n|/:e-*-^«(*-^)£e-*-^«^(l - 15)^(0)) 

(16) 



and 



CA^^[V,Al LoA=^[Ho,A]. (17) 



The choice of the initial conditions Eqs.(|6]),(l7]), and 
Eqs.jni) yield 



(18) 



N 



J2Ut) = l['\t) = - ^ /(2)(t). (19) 



n=2 



(2) 

Evidently, the definition of /„ (t) is crucial for the calcu- 
lation of the GME dSJ . For li^^ (t) we obtain the following 
expression (see details in Appendix) 

I^\t) = Re[ dTWin{t~T)B, n = 2,....Ar(20) 

JO 

2i SfoM^ 



B = exp 



a—h,v 



- 1 



where the variable S"q (see Eq. (jA10[) ) is determined by 
the reorganization energy of the a vibrational mode E^q, 
Eq. (fT2|) . when electron moves from the ground state |0) to 
the excited state |1) of the dimer P. As a check, using the 

(2) 

same techniques, we have calculated (t) and obtained 
the fulfillment of Eq.p9l). 

Note that if one considers the contribution of the non- 
diagonal density matrix elements, the first order term 
contributes to the GME as weU 123. [23. 



III. MODEL OF PRIMARY STAGE OF 
ELECTRON TRANSFER IN A REACTION 
CENTER 

To demonstrate the viability and utility of our ap- 
proach, we consider the electron transfer in RCs within 
a three-site model. In the RC, after photon absorption 
at the bacteriochlorophyll dimer molecules (molecule 1) 
the electron transfer may occur either through the M- 
side bacteriochlorophyll (molecule 2) or the L-side bac- 
teriochlorophyll (molecule 3) accessors. We assume that 
there are two non-zero coupling terms V12 and V13 and 
a forbidden direct electron transfer between molecules 
2 and 3 (V23 = 0). The two possible electron trans- 
fer pathways are related by the C2 symmetry axis. The 
pathway symmetry is broken by differences in amino acid 
around the donor-acceptor pairs in the different branches. 
These differences inhibit a charge separation through the 
M-branch in the wild type RC. The preference of the 
pathway (1 3) relative to the pathway (1 — >■ 2) is as- 
sumed to be governed by the energy differences and by 
the different couplings between donor-acceptor pairs in 
the branches. Indeed, experimental and theoretical es- 
timates for the energy differences provide the following 
figures [13, [H]: £1-63 « 0.05 -0.1, £2-^1 « 0.1 -0.2 
(eV). 

The Hamiltonian Hq determines the vibrational state, 
the equilibrium position of the oscillatory mode, and, 
consequently, the density operator for the vibrational 
subsystem, when an electron is in the base state of the 
donor molecule before the excitation. The phonon bath 
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FIG. 1: (Color online) Time dependence of the site occupation 
probabilities Pi{t). Solid (dashed) lines are used for results 
with (without) the ITs. The basic parameters are used for the 
both panels. The specific parameters are: Sio = l,Sio = 0; 
a)2r2/fi = 2T3/h = (5ps)-^ b)2r2//i = 2r3/ft = (Ips)-i 



is described by two vibrational modes, high frequency 
mode ujfi and low frequency ujy. The low- frequency mode 
characterizes the exterior medium phonon mode and the 
high-frequency mode describes the molecular vibrational 
modes of the donor and the acceptor centers [26|, |27|] . We 
take into account a relaxation time for the vibrational 
mode i /tq with the aid of a phenomenologi- 

cal parameter Tq (see also [28]). As a result, for numerical 
analysis the variables (jAll[) . (|A12[) are considered in the 
form 

JVa = (2nQ + 1)[1 - e"*/^"" cos{uJat)] + ie"*/^"" sin(wa<) 

Ada = <3xp{~t / Tpa) sm{u!at) ~ exp{~T / Tpa) sin(a;Qr) 

To illustrate the contribution of the ITs in the ET, we 
use parameters that may elucidate in the Markovian ap- 
proximation the observed L-side experimental kinetics of 
wild-type (WT) RCs of sphaeroides [Hi^. Note 
that the probability of the M-side electron transfer was 
excluded in such considerations. 

In partiuclar, the following set of parameters (set I) is 
used to characterize the electron transfer in a wild type 
RC via the L branch (the pathway 1 — > 3): V13 = 7 meV, 
huJh = 187 meV, hiVy = 12.5 meV, S^'^ = 0.5, S^^ = 8 
[26(1 . The same set can be used for the electron transfer 
via the M branch (the pathway 1^2), taking into ac- 
count that there is a two-fold symmetry. To provide the 
asymmetry in the electron transfer via the two pathways 
we consider: — ei = —56 meV and £1 = 110 meV 
(see discussion in [13] )■ This choice forms the set II. 

Evidently, that the RC is an open system which in- 
teracts with another part of overall system. This part 
can be assumed to have a quasi-continuum spectrum. To 



mimic the realistic situation we introduce sink parame- 
ters Ti/h (i = 1, 3) which characterize the electron trans- 
fer to another part of the overall system with a quasi- 
continuum spectrum. This is an effective approach to de- 
scribe resonance scattering phenomena in open systems 
with a weak coupling to the environment (see, for exam- 
ple, [13]). When an electron is transferred to this sub- 
system, the backward electron transfer can be neglected. 
The primary charge separation step (1 — 3) occurs in 
purple bacterial RCs with a lifetime ~ 2ps at room tem- 
perature. The next step of electron transfer occurs in 
the time scale Ips [31|. In the numerical calculations 
we have used the following values of sink parameters: 
2ri/;i = (250ps)-i|3l, 2r2/;i = 2r3/fi = (5ps)-i and 
{lps)~^. The parameter 2Ti/h characterizes the decay 
of the system to the ground state, while 2T2.3/h are as- 
sociated with the electron transfer to the next molecules 
(subsystems) which are beyond the scope of the present 
analysis. The flow direction also depends on the param- 
eter which characterizes the amount of energy stored 
in unrelaxed high frequency vibrational modes. The scale 
reorganization constant values are chosen as Siq = 0.5 
and = 1. We take Ta — 3.5ps for the lattice re- 
laxation time, in accordance with the observation that a 
vibrational mode relaxation time is of order a few ps [s^ . 
The set III consists of 2ri/;i and t. The sets I, II, III 
form the basic parameters of calculations of occupation 
probabilities at 300K, shown at Figs. 1,2. 

Let us consider a case when the lattice relaxation time 
is smaller than ones defined by the sink parameters (see 
Fig. la). The system has time to achieve the regime when 
the impact of the unrelaxed phonon mode on the occu- 
pation probabilities is easing. This imposes the impact 
of the ITs upon the electron flow direction. The electron 
flow direction is similar to the one defined by the theory 
without the ITs. At the early stage of electron transfer 
the ITs have, however, a strong influence on the quan- 
tum yields of the electron flow through different branches. 
The theory without the ITs predicts that the M branch is 
inactive. The ITs activates the electron escape through 
the branch M (molecule 2). 

The energy stored in the unrelaxed phonon modes is 
large {Siq — 1) in this case. The parameters h/2V2 = 
h/2T^ = bps are large enough in comparison with the 
phonon relaxation time Tp = 3.5ps. The system evolves 
in time to the regime where the ITs influence decreases on 
the electron transfer (the system begins to forget the ini- 
tial conditions). Therefore, the flow direction is the same 
as in the case without the ITs. If the sink parameters are 
smaller than the phonon relaxation time (see Fig. lb) the 
ITs change the favourite partway for the electron trans- 
fer: the dominance of the pathway (1 — 3) (without the 
ITs) is replaced by the dominance of the pathway (1 — 2) 
(with the ITs). 

Let us consider the regime, when the unrelaxed phonon 
modes of the medium are taken into account (6*^0 2) 
(Fig. 2). Such amount of the unrelaxed medium modes 
has no a strong impact on the main characteristic of elec- 
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FIG. 2: (Color online) The basic parameters are the same as 
in Fig. lb. The specific parameters for this figure are: Siq = 2. 
a)S'i'o = 1; b)5'?o = 0.5. 



tron transfer. The medium modes can not store large 
amount of the unrelaxed energy. Indeed, the transition 
of the system to the excited state affects rather the rear- 
rangement of atom positions in the donor-molecule and, 
thus, has no a strong impact on the medium atom posi- 
tions. The ITs increase the probability for electron trans- 
fer via the molecule 2 in comparison to the case with- 
out the ITs. The electron lifetime in the system is still 
short enough not to forget the initial condition. With 
the decrease of unrelaxed phonon energy in the system 
(compare Figs. 2a, b) the importance of the ITs is also 
decreasing. 

Another important ingredient is the dependence of the 
results on the ratio of the phonon relaxation time and 
the time associated with the interaction V. All results 
discussed above are related to the situation when r > 
fi/Vim- In this case the electron transitions are fast, while 
the phonon relaxation is a slow process. The ITs will 
contribute to the ET depending on the amount of energy 
stored in the unrelaxed vibrational modes, in accordance 
with the discussion above. On the other hand, if r <C 
fi/yim the ITs produce a marginal effect, and the initial 
conditions are forgotten relatively quickly. 



IV. SUMMARY 

In conclusion, we suggest the microscopic approach to 
study the effect of the ITs on the electron transfer in a 
system with initial conditions being far from equilibrium. 
The IT impact depends mainly on the amount of energy 
stored in the initially unrelaxed phonon modes and also 
on the lifetime of the electron in the system. If the elec- 
tron lifetime is much longer than the phonon relaxation 
time, the ITs do not affect the quantum yields of elec- 
tron transfer via possible pathways. In systems, where 



this condition is not fulfilled, the ITs can cause the elec- 
tron transfer via channels which are closed in the case 
without the ITs. 



Appendix A: Initial condition terms 

Let us employ the initial conditions ([51), ©, in the def- 
inition of the second order of the initial term p3|) .([TC )) . 



li^\t) ^ J* dTTr(^\n){n\Ce~'^°(''^^ X (Al) 
X £e-^^«^(pi-po)|l>(l| ] , n^2,...N. 



With the aid of Egs. pT)) we obtain 

12 rt 



dr X 



(A2) 



X exp ( --^i^n - ei){t - r) ) X 



X exp(^--(r„+ri)(i-r)j X 

X Trg fe*^i*e-^^"(*-^)e-*^i^(pi - 



Po, 



To proceed further it is necessary to calculate the trace 
over the environmental (vibrational) degrees of freedom, 
which has, in general, the following form 



I = Ttq (e^^'"*e-^^"(*~^)e-^^">fc) (A3) 
In the path integral formalism this trace can be written 

as 

oo noc 



I = / dq dx 



oo — oo 

oo poo 



oo 'J —oo 



X dy dyi{kva\q) 



(A4) 



X (g|exp( |xi)(a;i|exp (^^i^(t-Tyj \y) 

X (y|exp( -*^r] \y,){y,\e^p{~l3H]l)\ki.^)/Z. 



Here, Z — Trq [exp(— /3 H^)] and a vibrational state is 
{q\kvc) = {z)^N X H^^ (z) exp {-z^/2) 



1/2 



, z = (g - dka)/io , 



where £o = ^/h/muJa is the oscillator length, and H^^^ is 
Hermitian polynomial of order Va- 
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In order to obtain an analytical result we suggest to we obtain 
calculate this expression in the integral path formalism 
[3^ . As a result, Eq. (|A4|) takes the following form 



dq 



dxi 



dy 



dyiKmiq, -t;xi) X 
(A5) 



X - exp -5,7,„A/; exp 2i5;j,feA^„ , (A9) 



where 



K,{x,t]y) EE {x\exp{~iH^t)\y) 



1 



— . , -B-^ exp < I „ 

w 2i7r£Q sincjai 2£gsina;c(i 

Jb = coswQt[(a; - d^a)^ + (v - diaf 
-2{x - dia){y - dia) . 



(A6) 
(A7) 



In virtue of the equation 
^Xm(x - do)cxp[-l3{m + l/2)huJa]Xm{y - do) = 

m 



1 



■ exp 



2£l sinh(;3fta;„) 



Y 27r£2 sinh(/3ncjc<) 

Ti = cosh(/3Sa;„)[(a; - dof + (y - d^f 
-2{x - do){y - do) 



(A8) 



A/; 



Knn/f>^a , (AlO) 

(2n„ + l)(l-cosw„(t-r))+ (All) 
i s'mujait ~ t) I 

sinujat ~ sina;„r . (-^-12) 



Taking into account the definitions (IA2[) . (fTTj) . with the 
aid of Eq. (|A9p , for m = 1 we obtain the result (^0)1 . 
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